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ABSTRACT 



Aims. The effects of atomic diffusion on internal and surface abundances of A and F pre-main-sequence stars with mass loss are 
studied in order to determine at what age the effects materialize, as well as to further understand the processes at play in HAeBe and 
young ApBp stars. 

Methods. Self-consistent stellar evolution models of 1.5 to 2.8 M Q with atomic diffusion (including radiative accelerations) for all 
species within the OPAL opacity database were computed and compared to observations of HAeBe stars. 

Results. Atomic diffusion in the presence of weak mass loss can explain the observed abundance anomalies of pre-main-sequence 
stars, as well as the presence of binary systems with metal rich primaries and chemically normal secondaries such as V380 Ori and 
HD72106. This is in contrast to turbulence models which do not allow for abundance anomalies to develop on the pre-main-sequence. 
The age at which anomalies can appear depends on stellar mass. 

Conclusions. For A and F stars, the effects of atomic diffusion can modify both the internal and surface abundances before the onset 
of the main-sequence. The appearance of important surface abundance anomalies on the pre-main-sequence does not require mass 
loss, though the mass loss rate affects their amplitude. Observational tests are suggested to decipher the effects of mass loss from 
those of turbulent mixing. If abundance anomalies are confirmed in pre-main-sequence stars they would severely limit the role of 
turbulence in these stars. 
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1 . Astrophysical context 

The recent availability of magnetic field data from instruments 
such as ESPaDOnS at the Canada-France-Hawaii Telescope 
and Narval at the Bernard Lyot Telescope has allowed stel- 
lar physicists to explore long standing questions in stellar 
physics. A flurry of recent studies have specifically focused 
on observing globally organized magnetic fields in intermedi- 
ate mass pre-main-sequence (PMS) stars in order to determine 



stars (|Donati et al. 


1997r.lAlecian et al.Ll2005: Wade et al 


2005; 


Catala et all I2007: 


Alecian et al., 2008a b; Folsom et al., 


2008; 


Alecian et all 2009 


). This is particularly interesting since A and 



B type stars have, at most, a very thin convective envelope, 
which is probably not sufficient to generate solar type, dynamo 
driven, magnetic fields. With this in mind, the favored hypoth- 
esis used to explain the presence of important magnetic fields 
(from 100G to lOkG) in A and B main-sequence (MS) stars is 
the fossil field theory. In this context, the magnetic field would 
either have originated in the molecular cloud from which the star 
was born, or would have been generated by a dynamo process in 
the core during the star's earliest evolutionary stages. To test this 
hypothesis, the above mentioned studies have concentrated their 
efforts on char aterizing magnetic fields in Herbig (HAeBe) stars 
dHerbid. [T960h. which ar e belie ved to be the PMS progenitors of 
Ap/Bp stars dWade et all[2005l) . 



In some cases, namely for HD72106, RS Cha and V380 Ori, 
binarity allows us to probe the effects of chemical separation 
since both stars likely had the same initial composition. In fact, 
for at least two of these systems, the secondary has a solar metal- 
licity, whereas the heavier primary star has an above solar metal- 
licity and is likely chemically anomalous. Since atomic diffu- 
sion timescales generally decrease when stellar mass increases, 
chemical separation could perhaps offer an explanation. 

Furthermore, a more thorough analysis may allow to charac- 
terize other phenomena at play within these stars. What can the 
observed abundances tell us on the various processes, such as 
convection, mass loss and magnetic fields, which compete with 
atomic diffusion? And on what timescales can we expect signif- 
icant surface abundance anomalies? This paper will adress both 
of these questions. 

InlVkk et alj d2010h (hereafter Paper I), stellar evolution 
models with mass loss were introduced and shown to reproduce 
observed surface abundance anomalies for many AmFm stars. 
However, observed abundance anomalies do not allow to deter- 
mine whether it is mass loss or turbulence which is competing 
with atomic diffusion within the radiative zone of these stars. It 
was nonetheless established that the surface anomalies of AmFm 
stars were modulated by chemical separation which occurs deep 
within the star. Indeed, for both the mass loss models and models 
with turbulent mixing dRicher et all l2000b [Richard etldl l200lL 
and references therein), chemical separation near log AM/M« 
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Fig. 1. The position in the HR diagram (a) is shown for four models with a mass loss rate of 5 x 10~ 14 M o yr _I . The curves end 
(on the left) at 30 Myr, age at which all the models are on the main-sequence. The dotted segment of the 2.80 Mq curve represents 
a model with a mass loss rate of 10~ 13 Moyr _1 and was added in (a) and (c) in order to facilitate extrapolation, though it was not 
added in other pane ls for which the ma ss loss rate has a n effect. Observatio ns are shown i n the HR diagram f or three sets of binary 
stars: (■) V380 Ori (lAlecian et al.ll2009t) : (a) HD72106 dFolsom et alJl2008h : (•) RS Cha (lAlecian et alJl2005h . The evolution of the 
bottom of the surface convection zone (b), T e g (c), as well as the abundances of He (d), 7 Li (e), Ca (f) and Fe (g) are also shown. In 
panel (c), (x) marks the approximate end of the PMS. 



-5 to -6 is responsible for the anomalous surface behavior. For 
PMS stars on the other hand, chemical separation which occurs 
at this depth cannot explain observed abundance anomalies since 
the timescales are much too long. In mass loss models, chemical 
separation occuring near the surface leads to anomalies which 



appear at the surface within a few Myr, and could therefore rec- 
oncile observations. 

In the following analysis, and in our calculations, mass loss 
is considered in non rotating stars, since in such stars, mass 
loss could be the only process competing with atomic diffusion 
within the stable radiative zones. An outline of the main aspects 
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of our evolution code is found in Sect.[2] The stellar models will 
be presented in Sect. [3] with particular attention to the various ef- 
fects of mass loss and atomic diffusion in the interior (Sect. 13. Il l 
and at the surface (Sect. 13. 21 of PMS stars. In Sect.|H the models 
are compared to observations of various HAeBe stars. Finally, in 
Sect.[5j the implications of our results on the processes involved 
in HAeBe and, by extension, in ApBp stars will be discussed. 

2. Calculations 

The detailed description of the stellar evolution code used for 
these computations can be found in Paper I and references 
therein. At age zero, all models are ful ly convective with th e 
abundance mix prescribed in Table 1 of iTurcotte et al.l (Q998). 
Opacities are continuously updated for every mesh point as 
abundances evolve. The only adjustable parameter, the mixing 
lenght parameter a, was calibrated by fitti ng the current radiu s 
and luminosity of the Sun (see model H of ITurcotte et aLp 998). 
Radiative accelerations are taken from Richer et al. (1998) with 
corrections due to re distribution from lGonzalez et all d 1995b and 
iLeBlanc et al.l (l2000l) . The introduction of mass loss and its im- 
pact on transport are discussed extensively in Sect. 4 of Paper I. 

No extra mixing is enforced outside of convection zones. The 
effects of atomic diffusion materialize as soon as radiative zones 
appear as the models evolve toward the MS. These effects be- 
come more important as the radiative zone expands toward the 
surface, where atomic diffusion timescales are much shorter. 

The unseparatecQ mass loss rates considered range from 
10~ 14 to 10~ 13 M Q yr~': mass loss rates which lead to surface 
abundances compatible with observations of many AmFm stars. 
Due to uncertainties related to the nature of winds for A and F 
stars (see discussion in Sect. 4.2 of Paper I), we have chosen to 
limit our investigation to unseparated winds in order to avoid 
introducing additional adjustable parameters. 

Finally, th i s pape r is part of a series of papers starting with 
ITurcotte et al. (1998), where the mixing length used was cali- 
brated using the Sun for given boundary conditions, helium and 
metal abundances. For consistency, and in order to isolate the ef- 
fects of the processes of interest, the same boundary condition, 
solar composition and mixing length are used for all calculations 
of Pop I stars (see also Sect. 2 of Paper I). 

3. Evolutionary models 

In Fig.Q] the position in the Hertzsprung-Russell (HR) diagram, 
the evolution of AMbsczB T e g as well as of the surface abun- 
dances of He, 7 Li, Ca and Fe are shown for four stellar models 
with the same mass loss rate (5 x 10~ 14 M o yr _1 ). Due to numeri- 
cal instabilities related to the complete disappearance of the sur- 
face convection zone, the2.8OW5E-140 model could not be con- 
verged up to the main-sequence. Therefore, in order to complete 
the HR diagram, the curve was continued using a 2.80W1E-13 
model. In this case, the effect of doubling the mass loss rate 
on the position in the HR diagram is smaller than the width of 
the line. Its effect on surface abundances will be discussed in 
Sect.ll2l 

The exact definition of the MS and, by e xtension, th e PMS 
is somewhat arbitrary. The definition given bv libenl(ll965l) stipu- 
lates that the zero age main-sequence (ZAMS) begins when the 

1 Having the same composition as the photosphere. 

2 BSCZ: bottom of surface convection zone. 

3 The expression 2.80W5E-14 corresponds to a 2.80 M G model with 
a mass loss rate of 5 x 1CT 14 M G yr~' . 



thermal-gravitational energy is reduced to one percent of the lu- 
minosity of the star, although this may d epend on such things as 
the treatment of convection dCoxl [19681) . According to this def- 
inition, immediately after of the end of the PMS, indicated by a 
(x) in Fig.QJ, T e ff and L begin varying much more slowly. The 
time spent on the PMS varies from ~4Myr for the 2.80W5E-14 
model to ~19Myr for the 1.50W5E-14 model. 

3.1. Radiative accelerations, internal abundances and 
structure 




-12 -10 -8 -6 -4 -2 



log tM/M, 

Fig. 2. Internal variation of V r - V ac | for three models with a 
mass loss rate of 5 X 10 _14 M Q yr _I at the onset of the ZAMS. 
The surface is at log AM/M„=-12 and transport is convective 
when V r - V a( j > 0. For different masses, the position of convec- 
tion zones in relation to log AM/M, changes slightly, though it 
is constant in relation to T. 



Even on the PMS, the internal structure and abundances vary 
significantly between 1.5 Mq and 2.8 M Q . In Fig. [2] convection 
zones are shown for models of different mass as they arrive on 
the ZAMS (see also Fig.Q}?). The internal distribution of con- 
vection zones is strongly correlated with stellar mass. For the 
1.5OM model, the SCZ includes the linked H and He con- 
vection zones and is never shallower than log AM/M* = -7.2 
(see also Fig.Q}?). For the 1.90W5E-14 model, the surface con- 
vection zones separate into a linked H-Hel CZ and a deeper 
He II CZ, while a separate He I CZ materializes in the heavier 
2.50W5E-14 model. All convection zones disappear completely 
in the 2.80W5E-14 model before it has even arrived on the MS 
(see Fig.Q}?). Finally, over the entire PMS phase, the opacity 

bump near log AM/M, 7.2 resulting from iron peak element 

accumulation is not large enough to induce an iron peak convec- 
tion zone for any of the three models, though iron peak opac- 
ity might extend the SCZ slightly inward for the 1.50W5E-14 
model. If the mass loss rate is < 10~ 14 M o yr _1 , an iron convec- 
tion zone may appear after the PMS (see Fig. 5 of Paper I). 

Significant abundance variations appear in the interior of 
our PMS models. For example, in the 2.50W5E-14 model (bot- 
tom row of Fig. [3}, a 1.1 dex Ca overabundance develops at 
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log AM/M. 

Fig. 3. Comparison of the normalized local flux with radiative accelerations and internal abundances of 7 Li, O, Ca and Fe at 3 
different ages for the 2.50W5E-14 model. The curves end (on the left) at the bottom of the surface convection zone. In the middle 
row, the dotted line represents gravity. 



log AM/M* ^ -9.4 as early as 4Myr. It might seem puzzling 
however that such an overabundance transforms into a strong 
underabundance which reaches below -1 dex at 30Myr. Why 
should the behavior at 4 Myr be so different from the rest of the 
star's evolution and why does a strong overabundance develop in 
a region where g ra d(Ca) is below gravity? The reason stems from 
the fact that the wind progressively advects matter from deeper 
within the star. That depth is simply given by: 

AM a Mt. (1) 

So, at 4 Myr, the wind brings to the surface matter which origi- 
nates from4xl0 6 yr-5x lO^Moyr -1 = 2xl0" 7 M o ~ lO^M*. 



Correspondingly, one sees in the top row of Fig.[3]that the flux is 
nearly constant from the surface down to that depth (except over 
CZs), while it is clearly not constant below that depth. Naturally, 
the depth above which the flux is conserved due to advection 
from the wind increases with age (compare the curves at 4, 30 
and 300 Myr). In order to conserve the flux, X{Ca) at 4 Myr in- 
creases above AM/M, ~ 10~ 7 to compensate for the decrease in 
grad(Caj3- This process can be described by: 

T(r) = est r 2 p(U + C/ W )c (2) 



4 The same type of solution was obtained for oxygen in 
lLandstreet et~alli 19981) . 



M. Vick, G. Michaud, J. Richer, and O. Richard: Abundance anomalies in pre-main- sequence stars 5 



i i i I i i i I i i i I i i i I I i i i | i i i | i i i | i i i I I i i i | i i i | i i i | i i i I I i i i | i i i | i i i | i i i I I i i i | i i i | i i i | i i i I I i i i | i i i | i i i | i i i 




log AM/M, 



Fig. 4. Wind velocities (dot-dashed line: 10~ 14 M o yr _1 ; dashed line: 5 x 10 _14 M Q yr _1 ) and diffusion velocities (solid when positive, 
toward the surface, and dotted when negative) of a few selected elements in a 1.90 M model near the ZAMS. For most species, 
wind velocities decrease more rapidly inwards than diffusion velocities. 



where T(r) is the local flux at radius r, p is the local density, 
U and £/ w are the advective part of the atomic diffusion velocity 
and wind velocity respectively, and c is concentratioifl Although 
gravity is stronger than g ra d over a fraction of the stellar envelope 
for all elements shown in Fig. [3] the downward diffusion velocity 
is never larger than the wind velocity (see Fig.|4]for an example 
with a 1.9 M model). Therefore, as long as the absolute value 
of the wind velocity is larger than the downward diffusion ve- 
locity, any given element is dragged toward the surface, and its 
local abundance adjusts in order to conserve flux. Notice that the 
fraction of the envelope AM which can be described by Eq. [[2) 
increases with time since the wind progressively advects more 
mass toward the surface. 

With this in mind, one can understand that at t — 4Myr, 
the increase of X(Ca) above AM/M* ~ 10~ 7 is generated by 
the large flux arriving from regions where g ra d(Ca) is larger than 
gravity and which is conserved even as g ra d(Ca) decreases to- 
ward the surface. As time passes, Ca arrives at the surface from 
deeper within the star, where g ra d(Ca) is much smaller than its 
value at AM/M* ~ 10~ 7 , so that the flux is smaller, the flux con- 
servation applies over a larger mass and the surface abundance 
consequently decreases. This shows the importance of a solu- 
tion over the whole star, as done here, since applying an inner 
boundary condition at AM/M, ~ 10~ 7 would lead to an erro- 
neous solution after 4 Myr. 

On the PMS, the internal concentration variations are much 
smaller for 7 Li, O and Fe than for Ca. However, after the star 
arrives on the MS, as illustrated by the curve at 30 Myr, larger 
variations appear, including a nearly ~0.8 dex overabundance of 
Fe which spans from the surface down to log AM/M, ~ -8. 

As shown in Fig.|4l for a mass loss rate of 10~ 14 M o yr _1 , 
the downward diffusion velocity is greater than the wind veloc- 
ity for some elements, whereby the flux conservation regime 
as described by Eq. [[2) cannot be extended to all elements. 
The regime shift approximately occurs at 2 x 10~ 14 M yr _1 (see 
Fig. 5 of Paper I and corresponding discussion). This mass loss 
rate also marks the limit below which iron accumulation near 
T ~ 200 000 K can lead to iron peak convection (see Fig. 5 of 
Paper I). 

3.2. Surface abundances 

An element's surface abundance depends on age, stellar mass 
and mass loss rate. In Fig.Q] the surface abundances of He, 7 Li, 

5 This equation is derived and discussed in Sect. 5.1.1 of Paper I. 



Ca and Fe are shown for models of 1.5, 1.9, 2.5 and 2.8 M 
with a mass loss rate of 5 x lO _14 M yr _1 . As stellar mass in- 
creases, anomalies appear at the surface earlier. For instance, at 
around 3 Myr the 2.80W5E-14 model has a 0.6 dex overabun- 
dance of Ca while all other models still have their initial abun- 
dances (Fig.QJ). The same can be said for Fe overabundances or 
He underabundances, which appear later in smaller stellar mass 
models. Nonetheless, it is evident that for a given element, the 
overall shape of the surface abundance evolution curve is very 
similar for the three heavier models. This is due to two things: 
the appearance of a separate He II CZ, as well as flux conserva- 
tion as described in the previous section. 

For all four elements shown in Fig.[T| the initial, short-lived 
abundance maxima are caused by evolutionary effects. The ap- 
pearance of a radiative zone which separates the He II CZ from 
the SCZ allows for chemical separation to occur near the sur- 
face. The direction of the anomaly is determined by g m d-g within 
that region (between log AM/M, ~ -10.5 and -9), and thus 
He, 7 Li and Ca become underabundant, whereas Fe develops 
an overabundance (compare middle panel of Fig.[3]with Fig.[TJ. 
Following this brief episode, the surface abundance is deter- 
mined by flux conservation, and so results from chemical sep- 
aration occuring deeper within the stafl For instance, the subse- 
quent, gradual X(Fe) increase results from the wind slowly ad- 
verting matter which originated below log AM/M, 7 (for the 

2.50W5E-14 model), where g ra d(Fe) is greater than g. Similarly, 
the variations of X( 7 Li) which materialize at the surface are due 
to variations seen in g ra d( 7 Li) between log AM/M, ~ -7.5 and 
-6. 

The amplitude of the anomalies also depends on stellar mass. 
The Ca overabundance that materializes on the PMS in the three 
heavier models reaches 0.65 dex for the 1.90W5E-14 model, 
while it reaches 1.1 dex in the 2.50W5E-14 model. The He and 
7 Li underabundances, as well as the initial Fe overabundance fol- 
low the same dependence. The 1.50W5E-14 model behaves dif- 
ferently since most of the mass interval from log AM/M, ~ -7.5 
and -6 is mixed by convection (see Fig. [2]). 



6 Thermohaline convection, as described in Thead o~et all j2009t> . 
should not affect our results as the small mean molecular weight gra- 
dients that materialize occur in or just below the He II CZ, and should 
not affect surface abundances — after a fraction of a Myr (see Sect. l3.lt 
— if the internal solution is dominated by the wind (e.g. M > 2 x 10 -14 , 
see also Fig. 9 and discussion in Sect. 8. 1 of Paper I). 



6 



M. Vick, G. Michaud, J. Richer, and O. Richard: Abundance anomalies in pre-main-sequence stars 





/ ~p -1 
X 

\ 

3 
m 

>< l 
o 





: 1 1 


4 He; 




1 - 
7 Li- 

^ : 




i - 

- 




1 i - 

s - 


cr- 








- ,f| 






— i — E 


E — i 1— 


/ g Mn- 


E — I 1— 

: A 


Fe- 




r«" : 




z - 


i 




■ 


1 




1 


; r 


i 



10 20 



10 20 



10 20 
Age (Myr) 



1 

1 


1 


Ca- 


1 


1 






lO" 12 
lO" 13 
5x 10" u 

10 -14 

I 




1 



10 20 



10 20 



Fig. 5. Surface abundance variations for models of 1 .9 M with different mass loss rates which are identified in the lower right panel 
in M yr~' . The vertical line indicates the approximate end of the PMS. 



3.2.1 . Effect of varying the mass loss rate 

Generally, increasing the mass loss rate lowers the amplitude 
of surface abundance anomaly extrema, as sinking elements are 
further advected toward the surface while supported elements 
are more effectively evacuated. In Fig. [5] the surface abundances 
are shown for a 1 .9 M model with four different mass loss rates. 
For this model, the PMS ends near lOMyr (see Fig.QJ). For the 
1.90W5E-14 and 1.90W1E-13 models, He, O and S decrease 
slowly and monotonically, while Li oscillates around its initial 
value. A calcium overabundance occurs at ~7 Myr, after which 
X(Ca) returns to its original value near the end of the PMS. Iron 
peak elements are overabundant for both models throughout the 
PMS and onto the MS, with Fe maxima of 0.4 and 0.6 dex for 
the 1.90W1E-13 and 1.90W5E-14 models respectively. For all 
elements, maximum amplitudes are smaller when the mass loss 
rate is larger. 

The 1.90W1E-14 behaves quite differently: large anomalies 
appear very rapidly, especially for sinking elements; 7 Li and O 
underabundances reach —1.2 and -2 dex respectively, while Ca 
is more than 100 times underbundant as soon as chemical sep- 
aration reaches the surface. This is because, in contrast to the 
other two mass loss rates, the wind generated by a mass loss 
rate of lO~ 14 M yr~' is not strong enough to dominate inward 
diffusion for most sinking elements. In Fig. [4] for instance, one 
can see that for this mass loss rate, the inward diffusion velocity 
of 7 Li, O and Ca is greater or equal to the wind velocity over 
an important fraction of the upper layers. On the other hand, 
supported elements behave as in the other two models, though 
with larger overabundances (at least up to the age of the last 
converged model). This is logical since the weaker wind cannot 
evacuate as much of the metal rich material which has accumu- 
lated in the SCZ. 

On the PMS, and for this stellar mass, only the 1.90W1E-14 
model generates sufficient anomalies to cause Z to vary signifi- 
cantly from its original value. The 1.90W1E-12 model only al- 
lows very small abundance anomalies at the surface (~ 0.1 dex). 

3.2.2. Comparison to turbulence models 

In the models of iRicher et al.l d2000l) . turbulent mixing down to 
T ~ 200 000 K, the depth of the iron peak CZ, is assumed to be 



the process which competes with atomic diffusion in the stellar 
interior. However, due to efficient mixing, abundance anomalies 
do not appear at the surface of these models until they have ar- 
rived on the MS. The signifi cant overabund ance of Ca in mass 
loss models (Fig. [If, see also lAlecianl fl996i. is not obtained in 
the above mentioned models with turbulence, and so could be 
used as an observational test for young stars. 

4. Comparison to observations 

Confronting our models to observations can offer constraints on 
mass loss rates as well as on mixing processes near the surface, 
such as convection. However, in order to compare models with 
observations, one must try to fit various parameters simultane- 
ously, including age, mass (or T e ff) and initial surface compo- 
sition. To reduce the arbitrariness of the comparison, we chose 
three binary systems, since both components of such systems 
should have the same age and initial abundances. Unfortunately, 
determining this age and initial abundance mix is difficult. For 
all three systems, the age deduced by observers was determined 
using model isochrones in which many assumptions were nec- 
essary. For instance, in all three cases atomic diffusion was ne- 
glected, and large mass loss or accretion rates (~ lO _8 M yr _1 ) 
were assumed. For this reason, we believe that it is preferable 
to fit the position in the HR diagram from this study rather than 
use the age determined by other isochrones. As for initial abun- 
dances, two of the systems have a chemically "normal" sec- 
ondary which offers a glimpse into the initial metal content, even 
though "normal" abundances can mean different things depend- 
ing on the solar mix used as a reference^. However, assuming 
that both components of a given binary system have the same 
initial abundances allows us to make a differential comparison 
of surface abundances and conclude that differences between 
the two are caused by internal processes such as atomic diffu- 
sion. The chosen original abundances aren't crucial for our anal- 
ysis since our objective here is to show when abundance anoma- 
lies appear on the PMS, which is not very sensitive to the ini- 

There is still much debate on the newer lAsplund et all J2005T) 
and lAsplund et al.l d2009h initial solar abundance mix for which abun- 
dances of CNO are significantly lower than the older solar mix of 
iGrevesse et aT Hl996l); the latter are in bett er agreement with helioseis- 
mology models TDelahave & Pinsonnea ult 2006). 
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tial abundances used to construct the models (see Fig. 13 and 
Sect. 8.1 of Paper I). 
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Fig. 6. Surface abundances for models representing the compo- 
nents of the binary star V380 Ori. The models were chosen at 
the age which best fits the primary in the HR diagram (Fig.QJj). 
In both models, the mass loss rate is 5 x 10~ 14 M Q yr~'. 



In Fig.Q] models are compared to three observed binary 
systems: V380 Ori (■), HD72106 (A) and RS Cha (•), which 
are all believed to b e HAeBe stars. For the Ae star V380 Ori 
dAlecianet al.ll2009l) . the 2.87 M primary and 1.6M chemi- 
cally normal secondary are compared to the 2.80W5E-14 and 
1.50W5E-14 models respectively. Within the error bars, the fit 
in the HR diagram is very good and the c orresponding age of 
the primary is ~2.7Myr (the age given in lAlecian et alj 120091 
is 2+1 Myr). The 2.80W5E-14 model has developed overabun- 
dances of Fe and Ca which reach ~ 0.3 dex and ~0.6dex re- 
spectively. Though this model could only be converged un- 
til 2.8 Myr, it has already developed anomalies which explain 
t he hi gh metallicity ([M/H]=0.5) determined in lAlecian et alj 
(2009). As illustrated in Fig. [6] most metals heavier than Z > 15 
are overabundant by a factor of ~2-3, while CNO are barely un- 
derabundant. Lithium becomes -0.2 dex underabundant. At the 
best fitting age for the primary, the 1.50W5.14 model still has 
its initial abundances, which also agrees with observations. The 
potential effect of a magnetic field, which is observed in the pri- 
mary, will be discussed in Sect.|5] 

For HD 72 106, the 2.4 M primary and the 1 .9 M Q secondary 
are compared to the 2.50W5E-14 and 1.90W5E-14 models re- 
spectively. In the HR diagram, the fit is not perfect, especially 
for the primary, tho ugh this was also problematic in Fig. 4 of 
iFolsom et al.l 12008). In fact, these authors suggest that at a deter- 
mined age of 6-13 Myr (our best fit model is between 5-7 Myr), 
the primary is most likely on the ZAMS rather than the PMS. 
Nonetheless, they found that the primary was chemically anoma- 
lous with important overabundances of iron peak elements and 
a strong underabundace of He, while the secondary is almost 
solar. In the best fitting models for the primary (5-7 Myr), iron 



overabundances vary from about 0.1 to 0.3 dex, while Ca over- 
abundances vary from about 1 dex down to about 0.3 dex. For 
the same age interval, the He underabundance varies from -0. 1 
to -0.15 dex. These amplitudes are smaller than the observed 
values. This could be due to the mass loss rate being too large 
(see Fig. [5}. The pre sence of a magnetic field and of phase vari- 
ations (see Fig. 9 of iFolsom et al.ll2008h do not justify trying to 
achieve a better fit. A precise model would require 2 or 3 D cal- 
culations. Furthermore, some uncertainty remains on the exact 
amplitudes since all their abundances were determined simul- 
taneously by fitting the observed spectra for a single averaged 
phase. Over this same age interval, the secondary still has its 
initial abundances, which agrees with observations. 

For RS Cha, the 1 .89 M primary and the 1 .87 M secondary 
may be compared to the 1 .90W5E- 14 model. In the HR diagram, 
the slight discrepancy with the model can be explained by the 
slight difference in stellar mass. Within reasonable error bars, 
the model is either chemically normal (up to 6.7 Myr), or has 
developed an overabundance of Ca of about 0.4 dex accompa- 
nied by small underabundances of 7 Li and He of -0. 1 1 dex and 
-0.05 dex respectively. This model does not seem to explain the 
iron enrichment factor of 1.5 obtained by the authors. However, 
in contrast to the two previous binary systems, both these stars 
have the same composition; therefore, there is no difference in 
composition between the stars to explain, and so the initial abun- 
dances could be responsable for the abundance anomalies with 
respect to the solar compositi on. The X-ray emission observa- 
tions of lMamaiek et al.l dl999) suggest that accretion could play 
a role. A smaller mass loss rate would also lead to larger anoma- 
lies. 

Other young single HAeBe stars which are not shown in 
Fig.Q]may also be compared to our models. The star HD 104237 
has a mass of about 2.3 M (r e ff = 8000 K), a luminosity of 
about 1.42L (Bohmetal. 2004) and an approximate age of 

2 Myr (Ivan den Ancker et al.lll998l) , as we ll as a magnetic field 
of about 5 G dDonatiet al.L Il997l 120001) . It is an HAeBe star 
for which lAcke & Waelkensl (120041) found that Si, Cr and Fe 
abundances were solar, which agrees with our results since 
even the heavier 2.50W5E-14 model is roughly normal until 

3 Myr (see Fig.[B. Similarly, HD 190073, which has an age of 
1.2+0.6Myr and a mass of 2. 85+0.25 Mp (r eff = 9250 K), was 
also found to be roughly solar (lAcke & Waelkens1l2004l) . Within 
the timescales shown in Fig.Q] this is also compatible with our 
results. Finally, though it has just recently embarked on the MS, 
the young magnetic Bp cluster star N GC 2244-344 whic h has 
r eff a 15 000 K and an age of 2 Myr dBagnulo et al.ll2004l) can 
also be compared to our results. A rough extrapolation sug- 
gests that the observed ~1 dex overabundances of Si and Fe and 
~2dex overabundances of Ti and Cr, as well as the ~1 dex un- 
derabundance of He could simply be the result of chemical sep- 
aration which began on the PMS. 



5. General discussion and conclusions 

Though it is often assumed negligeable for A and F type PMS 
stars, chemical separation resulting from atomic diffusion can 
affect both the surface and interior of these young stars. The 
mass in convection zones (Pig.[2j and internal concentrations 
(Fig-Ell can be modified before stars arrive on the MS. The am- 
plitude of the internal concentration variations depends on stellar 
mass. Equivalently, the age at which abundance anomalies ap- 
pear at the surface also depends on stellar mass. In the presence 
of weak mass loss, and for models with no turbulent mixing, ro- 
tation or magnetic fields, significant internal variations and sur- 
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face anomalies (both exceeding factors of 10 for some elements) 
appear as early as ~2 Myr in early-type A stars, and ~20-25 Myr 
in cooler F stars. 

Mass loss only slightly affects the age at which abundance 
anomalies occur at the surface, although its inclusion in our cal- 
culations was necessary in order to evolve models to ages of 
interest. This being said, mass loss does affect the maximum 
amplitude of surface anomalies (see Fig. [5}; a mass loss rate 
> 10 -12 M Q yr -1 nearly flattens surface abundance anomalies, 
whereas mass loss rates < lO _I4 M yr _I allow strong under- 
abundances of elements which are not supported by the radiation 
field, and may allow the appearance of an iron peak CZ. 

The important Ca overabundance materializing in our mod- 
els of PMS stars, which is not obtained in similar models with 
turbulence, could help distinguish between these two scenarios 
(see Sect. 13. 2. Zt . If abundance anomalies observed at the surface 
of A and F type PMS stars are in fact due to chemical separa- 
tion, this strongly favors the mass loss model as presented in this 
paper and should be confirmed by asteroseismology (see also 
Sect. 6 of Paper I). More precise determinations of individual el- 
emental abundances may also allow further differentiation. 

Atomic diffusion coupled with a mass loss rate which is 
compatible with observations of AmFm stars may elucidate 
why there are PMS binary systems such as V380Ori and 
HD72106 for which the primary is chemically anomalous while 
the secondary remains roughly normal (see Sect.|4]i. Since dif- 
fusion timescales vary more rapidly with stellar mass than with 
mass loss rate (compare Figs. Q] and |5j, smaller mass loss rates 
(or atomic diffusion on its own) may also explain these sys- 
tems. Observed abundance anomalies in other single PMS stars 
and young ApBp stars may also result from atomic diffusion 
(Sect.|4]i, without the need for other more exotic explanations. 
This being said, many PMS stars have observed magnetic fields, 
strong accretion rates (from stellar disks) or both. These phe- 
nomena could also have an effect on transport in the external 
regions of the star. 

A magnetic field can impact chemical transport through its 
effects on rotation via magnetic braking, by modifying convec- 
tion in the atmosphere (Catt aneo et al.ll2003|). as well a s by af- 
fecting atomic diffusion velocities (lAlecian & Stiftl2006l) . In sin- 
gle stars for which rotation cannot be slowed by tidal forces, 
strong magnetic fields may be the only process which suffi- 
ciently reduces rotation to allow chemical separation. Although 
it has to our knowledge never been shown that magnetic fields 
can completely eliminate convection, their topologies can in- 
duce horizontally depend ent convection and abundance profiles 
(Babel & Mich audlll991l), in addition to anisotropic mass loss 
and/or accretion (Thea do et al.l 12005). If mass loss is in fact 
anisotropic, t hen strong mass loss ra tes, as suggested for young 
PMS stars bv lBoehm & C atala 1995, may be present only where 
field lines are vertical, and thus be compatible with anisotropic 
surface anomalies. Similarly, magnetic fields may extend the up- 
per limit of y eq = 100 km s -1 , which was found to eliminat e 
the effects of atomic diffusion dCharbonneau & Michaudll 19881) . 
This may allow ano malies to develop even in the fast rotators 
among HAeBe stars (iDavis et al.lll983l) . 

However, until these processes are better understood, and 
magnetic field geometries are better constrained for individ- 
ual stars, it would be ill-advised to introduce additional mag- 
netic field related parameters into our calculations. The effects 
of magnetic fields on chemical transport may best be grasped 
through comparative mapping of superfi cial magnetic fields 
and abundances as done for roAp stars in iKochukhov & Wada 



(l2010h andlLuftinger et al.1 (l2010h . as well as for HAeBe stars in 
iFolsom et al.l (120081) . 

Accretion also affects elemental distribution in the atmo- 
spheres of stars. For accretion rates greater than approximately 
10" 12 MQyr _1 , the abundance profiles in the atmosphere, and 
in the interior, simply reflect those of the accreted material 
dTurcotte & Ch arbonnea uTl993l) . For accretion rates a few times 
larger than lO _14 M yr _1 , the atmospheric abundances can pre- 
serve the accreted material signature while allowing for abun- 
dance gradients due to chemical separation to develop below the 
SCZ. For even smaller rates, atmospheric abundances essentially 
reflect the result of chemical separation as if there were no ac- 
cretion. In any case, as soon as the star stops accreting, chemi- 
cal separation resulting from atomic diffusion dominates within 
1 Myr. 

In terms of stellar modelling, HAeBe stars are extremely 
complex: stellar winds, accretion, rotation and magnetic fields 
complicate simulations. Current models cannot account for all 
these processes without invoking multiple parameters which 
may blur any actual physics taking place within these stars. 
Though there are other processes involved, it is shown that 
atomic diffusion can lead to abundance anomalies on the PMS, 
and that neglecting its effects could have an impact on calibrat- 
ing atmosphere models. 
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